Ferroelectricity in strained Cao.sSro.sTiOa from first principles 
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We present a density functional theory investigation of strained Cao.sSro.sTiOa. 

We have determined the structure and polarization for a number of arrangements of 

Ca and Sr in a 2x2x2 supercell. The a and b lattice vectors are strained to match 

the lattice constants of the rotated Si(OOl) face. To set the context for the CSTO 

study, we also include simulations of the Si(OOl) constrained structures for CaTiC>3 
O 

and SrTiC>3. Our primary findings are that all Cao.sSro.sTiOs structures examined 

j> ■ except one are ferroelectric, exhibiting polarizations ranging from 0.08 C/m 2 for the 

o : 

lowest energy configuration to about 0.26 C/m 2 for the higher energy configurations. 
We find that the configurations with larger polarizations have lower c/a ratios. The 
net polarization of the cell is the result of Ti-0 ferroelectric displacements regulated 



o 
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by A-site cations. 



There is great interest in combining ferroelectrics with semiconductors. Potential de- 
vices include nonvolatile memory, reprogrammable logic, and even quantum computation^. 
Ideally, ferroelectrics would be grown on the technologically dominant Si(OOl) surface, and 
there has been great effort at growing SrTiOs on Si(001)^~-. The in-plane lattice parameters 
of SrTiC>3 are 1.7% larger than those of the (110) and (110) directions of the Si(001) surface, 
and thin films of compressively strained STO have been shown to exhibit the desired ferro- 
electric behavior-^. However, as the thickness of the film increases, the STO relaxes to the 
unstrained state and the ferroelectric behavior is no longer observed^. One possibility to 
induce ferroelectricity into SrTiOs films is to introduce dopants such as either Ba or Ca£~-. 
For example Ca(i_ x )Sr x Ti03 offers the possibility of invoking ferroelectricity in SrTiOs while 
at the same time maintaining a closer lattice match to the Si(001) surface. Although there 
have been a number of experimental studies aimed at identifying the crystal structure and 
lattice properties of unstrained Carx-^Sr^TiO^"— , to date the authors are not aware of 
any theoretical studies of the strained case. 

In this work we present a theoretical study of ferroelectrically induced polarization in 
bulk, strained Cao.sSro.sTiOs. Structural optimizations were done using the VASP code^^ 
with the Projector Augmented Wave method to treat the electronic structure problem^. 
The gradient corrected XC-functionals were treated within the Perdew, Becke and Ernzerhof 
scheme^. In all cases we used a supercell consisting of 2 x 2 x 2 primitive cells or two layers 
of four primitive cells each. This requires eight alkaline-earth atoms in the unit cell while 
allowing for the inclusion of rotation and tilting effects in the oxygen octahedra. Periodic 
boundary conditions were employed to simulate a spatially extended material. We used a 
4x4x4 k-point mesh and a planewave cutoff energy of 350 eV. The unit cell was strained 
to match the lattice constants of the 45° rotated Si(001) face. For our unit cell the a and b 
lattice vectors are 7.728 A. The c lattice vector is obtained by calculating the total energy 
for a range of values separated by 0.01 A approximately centered on the cubic value until a 
minimum energy configuration was bracketed. For the total energy calculations the lattice 
parameters were frozen but the atomic positions were free to relax until the forces were 
smaller than 0.005 eV/A along any cartesian direction for any atom. 

For the calculation of polarization we used the Abinit code^ and FHI98 pseudopotentials 
of the Trouiller- Mart ins type 18 . The structure was imported directly from the VASP code 
structural optimization step and no relaxation was performed for the lattice or ions. 
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FIG. 1. (Color online) Top view schematic of the configurations considered in this work. The left 
column is an ID for each configuration, the second column is a schematic of the first layer, the 
third column is a schematic of the second layer. The final column gives the relative energies. The 
white/red circles represent the positions of the Ca/Sr atoms on the xy-plane lattice face. 



In the chosen unit cell there are 8 A-sites (alkaline-earth) of which four are occupied 
by Ca and four by Sr. Most of the 70 possible arrangements of the cations are related by 
rotational and mirror symmetries. We examined the 9 configurations shown in Fig. [TJ 

Our primary result is that all the configurations we considered except one exhibit a net 
ferroelectric dipole moment. The lowest-energy configuration consists of alternating stripes 
of Ca and Sr running diagonally across the strained face of the perovskite. This configuration 
possesses a net polarization of 0.077 C/m 2 . We find that there is a trend for the lowest energy 
configurations to have the smallest polarizations and the largest c/a ratio's. We summarize 
this information in table HI 

Several studies on SrTi0 3 have shown that quantum fluctuations may compete with 
ferroelectric distortions, leaving the perovskite in an unpolarized state referred to as an 
incipient ferroelectric^ 1 ^. To address this possibility in the doped perovskite, we have 
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Configuration 


AE[eV] 


a ra ti° 


Polarization [C/m 2 ] 


Hollxyl2ol2xy 


0.000 


1.043 


0.077 


Hollxyl2yl2xy 


0.090 


1.048 


0.035 


llollxllxyl2xy 


0.094 


1.048 


0.037 


Hollxl2yl2xy 


0.304 


1.020 


0.000 
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0.309 


1.030 


0.250 
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1.030 


0.245 


Hollxyl2xl2y 


0.316 


1.033 


0.253 


llollxllyllxy 


0.331 


1.030 


0.268 


llollxllxyl2y 


0.385 


1.027 


0.031 



TABLE I. Tabulated data for the configurations at 50% Ca, 50% Sr pictured in fig. [TJ Included 
are the configuration identification, the total energy, relative energy and polarization for the lowest 
energy of each configuration. 



calculated the energy difference between the polarized and unpolarized cells for the two 
cases corresponding to the lowest energy configuration (Uollxyl2ol2xy) and the alternating 
stacked sequence of Ca and Sr (llollxllyllxy). We find the energy differences to be 0.07 
eV and 0.03 eV per formula unit (5 atom cell) for the lowest energy configuration and the 
stacked configuration respectively. If we use formula (1) from reference^ to estimate the 
contribution of quantum fluctuations we come up with about 0.003 eV/oxygen or about .01 
eV/formula unit. Thus it is likely that the polarization in these cells persists. 

We find that polarization tends to increase with decreasing c/a ratio. There is an excep- 
tion to this trend for the configuration labeled Uollxl2yl2xy which has no net polarization 
but has the smallest c/a ratio of those considered. We have plotted the relative energies of 
the configurations as a function of c/a ratio in figure HJ In this figure we label four fami- 
lies of curves which are differentiated based upon the AFD rotations and/or tiltings. The 
four distinct combinations of rotation and tilting are presented in figure E] which consist 
of projections down the (100), (010), and (001) faces of the Hollxyl2ol2xy, Hollxyl2yl2xy, 
Uollxllxyl2x, and Uollxl2yl2xy configurations respectively. For simplicity of expressing the 
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FIG. 2. (Color online) Plot of the relative energy as a function of c/a ratio for the configurations 
considered in this work. The curves were shifted so that the globally lowest energy configuration 
aligned with 0.00 eV. The labels I, II, III, and IV correspond to arrangements of a + a + c~, a + a + c°, 
a~b + c~, and a~b + c + tilting and rotation of the oxygen octahedra. 



rotations and tilting of the oxygen octahedra we employ the Glazer notation 21 . This notation 
was developed for ABX3 perovskites and assumes rigid rotations of the oxygen octahedral 
cages whereas we are considering the case of AA'TiOs and allow distortions of the oxygen 
cages. 

Analysis of the displacement patterns shows that the tetragonal strain is driven by the 
AFD octahedral rotations hybridized with the ferroelectric displacements (FE) rather than 
by the latter alone. Earlier theoretical work has shown that the introduction of dopants 
with smaller atomic radii into the A-sites is able to drive the tetragonal strain^ 2 .. Thus the 
distribution of Ca-cations determines the rotation of the oxygen cages versus the distortion 
along the tetragonal (c) axis. On the other hand, the polarization is clearly associated with 
the FE displacements. We find that FE displacements of the Ti atoms relative to the oxygen 
cages do not vary widely despite the range of polarizations for the different configurations. 
For example, the Ti atoms undergo an average FE displacement from their centrosymmetric 
positions in the oxygen octahedra of 0.213 A for the Uollxyl2ol2xy configuration. This 
configuration has a polarization of 0.077 C/m 2 , whereas the llollxllyllxy configuration Ti 
atoms have an average FE displacement of 0.134 A while exhibiting the one of the largest net 
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polarizations of 0.268 C/m 2 . Instead, we find that displacements of the A-site cations plays 
a significant role in determining the overall polarization of the cell. In the Uollxyl2ol2xy 
configuration, the A-site Ca cations have very large displacements along the c-axis. The 
displacement of these cations is such that the Ca-0 dipoles are anti-aligned with the Ti- 
associated polarization of the cell. If the Ca ions are shifted to be coplanar with the oxygen 
atoms in the same atomic layer thus removing their contribution to the polarization in 
the cell, the polarization increases from 0.077 C/m 2 to 0.186 C/m 2 . If all A-site cations are 
shifted to their nonpolar locations, the polarization increases further to 0.222 C/m 2 . In their 
optimal locations, the Ca-ions have large displacements of up to 0.589-0.650 A in the three 
lowest energy configurations. Other configurations have Ca-ion displacements of 0.000 to 
0.243 A. The resulting picture is such that the configurations with the lowest energy, groups 
I and II in Fig. [31 have the highest tetragonal strain and lowest polarization as featured by 
the relatively strong a + a + -type octahedral rotations. These rotations will be referred to as 
tiltings hereafter in the text. 

We consider the effect of the position within the cell of the atomic arrangements on the sta- 
bility and polarization by examining the three cases of the lowest energy Uollxyl2ol2xy con- 
figuration, the Hollxyl2xl2y configuration which may be thought of as a rock-salt type struc- 
ture, and finally the Uollxllyllxy configuration which is constructed of alternating planes 
of Ca and Sr ions in the (001) plane. Note that for the Hollxyl2ol2xy and Uollxyl2xl2y con- 
figurations, the sequence of the A and A' cations (A = Ca, A' = Sr) along the (a, b) lattice 
directions is A-A'- and so on, while that of the Uollxllyllxy is either A-A- or A'- A'-. Along 
the c lattice direction, the Uollxyl2ol2xy give either a stacking of A-A- or A'-A'- whereas 
the other two give A-A'- and so on. Considering these cations as hard spheres, the packing 
of A-A'- is more compact than that of A'-A'- in the ideal tetragonal structure, since the Ca 
ion is smaller than the Sr ion. We have confirmed that the Uollxyl2ol2xy and Uollxyl2xl2y 
undergo smaller compressive strains of (0.5%, 0.5%), and (0.9%, 0.7%) compared to (1.2%, 
0.7%) for the Uollxllyllxy configuration along the a and b lattice directions. This results 
in reduced polar AFD rotations. For the Uollxyl2ol2xy configuration, the average octa- 
hedral rotation about the tetragonal axis is 5.7°, while the Uollxllyllxy configuration has 
an average octahedral rotation of about 7.7°. In the c direction, both the Hollxyl2xl2y 
and Uollxllyllxy are more compact than the Uollxyl2ol2xy which restricts the range of 
displacements that the Ca ions are permitted. This in turn reduces their ability to respond 
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FIG. 3. Ball and Stick models of the lowest energy configuration corresponding to each combination 
of AFD rotations and tilting of the oxygen octahedra exhibited by configurations investigated in 
this work. We present views down the (100), (010), and (001) planes respectively. The models 
correspond to the configurations lying at the bottoms of the curves labeled in figure [TJ In order 
from smallest to largest balls, the O are red, Ca white Sr green, and Ti yellow. 



to the internal electric field established by FE distortions. 

To assist in understanding how doping affects the stability and polarization of the per- 
ovskite we have performed a series of calculations in which we include only the effects of 
FE displacements, FE displacements coupled with AFD rotations, and FE displacements 
coupled with octahedral tilting about one axis. The results are presented in table [Til We 
find that freezing out the tilting and allowing only AFD rotations of the octahedra coincides 
with an increase in the cell polarization relative to the case where only FE displacements 
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are allowed, contrary to what is found for STO^. This behavior is the result of large ~0.5 
A displacements of the Ca atoms out of the (001) plane of the oxygens in the case where the 
AFD distortions are suppressed. The Ca atoms are displaced such that the Ca associated 
electric dipole opposes the contribution from the Ti associated polarization. This is facili- 
tated by a large c/a ratio for the structure. When AFD rotations are included the strain on 
the octahedral cages is reduced and the tetragonal distortion is smaller. In this case the Ca 
atom displacements are reduced to about 0.1 A due to the steric hindrance of the oxygen 
atoms. Thus their associated electric dipole is of a smaller magnitude. 

The inclusion of tilting without AFD rotations also acts to reduce the cell polarization. 
Due to the compressive strain imposed in the plane, the oxygen octahedra are compressed 
and there is a large tetragonal distortion out of plane which allows for displacement of the Ca 
atoms. The distance between the Ca atoms and the nearest oxygen atoms in the Ti-0 plane 
above them (relative to oxygen octahedra that are displaced downwards due to polarization 
effects) is about 2.43 A as compared to the sum of the ionic radii of Ca and O which is about 
2.38-2.42 A assuming Shannon ionic radii 24 . Thus they are bounded in their displacement 
primarily by the size of the O and Ca ionic shells. For comparison, the sum of the ionic 
radii of Sr and O is 2.56-2.60 A. Generally speaking, the A-site cations exhibit distances to 
nearest neighbor oxygens from adjacent Ti-0 planes that correspond to within 0.05 A of the 
sum of their ionic radii. 

Based on these observations, the placement of the Ca and Sr atoms in the lattice can 
be explained. Fig. H] consists of two sets of schematics indicating only the directions of the 
displacements of the A-site planar oxygen atoms relative to the A-site cations due to tilting 
of the octahedral cages. Note that the four sets of curves in fig. [2] reduce to the two sets 
of oxygen displacement schematics (including rotation about the c-axis gives the complete 
4 sets of curves). A study of the figure indicates that once the displacement of the oxygens 
about one cation have been established, the displacements of all oxygens in the unit cell is 
specified, assuming rigid oxygen octahedra. 

In layer 1 (LI) of the fig. Hltop set of schematics, the schematic labeled "o" indicates that 
all of the nearest neighbor oxygens that are co-planar with the A-site cation are tilted towards 
the cation. This set of tiltings corresponds to the large out-of-plane displacements that have 
been observed for the Ca atoms in the simulations. The lowest energy configurations contain 
this tilting arrangement. This set of tiltings, coupled with the large A-site cation out-of-plane 



LI 



L2 



I/II 









/ 




y 


\ s 


o 




o 


\ 


o 


/ 


o 

Jf \ 



iii/iv ^y 



xy 



v y ' 


A 


xy 


\ 


\ 


y 


X 


xy 


\ v 


X 




y 






o 




o 






o 




o 


• v 


V 




• 











o 



xy 



V 






V 




V 


V 


o 




o 




o 




o 


• 






' 




• 


• 



O ' ' X ' o x r 

FIG. 4. Schematic diagrams of the (OOl)-plane projected tilting of the oxygens surrounding the A- 
site cations. The top set of eight panels show the direction of oxygen displacements corresponding 
to curves I and II of figure [2j The bottom set correspond to those of curves III and IV. In each 
case, the left 4 panels are for layer 1 (LI), the right set are for layer 2 (L2). Each schematic is 
labeled with an o, x, y, or xy to indicate its proximate placement in the unit cell. Cells adjacent 
to each other share oxygen atoms and direction vectors indicate that. 



displacements, is associated with the lowest energy configurations. In all 9 configurations 
studied here, no such set of tiltings is associated with the Sr atoms due to their large ionic 
radii. The tilting drives the A-site cations out of the plane. In order for the Sr atoms 
to displace, the tetragonal distortion must increase with a corresponding increase in strain 
energy. On the other hand the Ca atoms are able to displace with no additional contribution 
to the tetragonal strain energy due to their smaller atomic radii. To quantify this, note that, 
in the Hollxyl2ol2xy configuration, the 0-0 nearest neighbor separation in A-site plane is 
3.29 A. The in-plane distance from one of the oxygens to the center of the square face where 
a non-displaced A-site cation would reside is 2.33 A. In order to accommodate a Ca ion, 
given the ionic radii of the Ca and O, the Ca ion would need to be displaced by about 0.59 A 
which is the same as the observed Ca displacement. In order to accommodate an Sr cation, 



10 

the tetragonal distortion would need to be increased by 0.53 A per layer. 

Given the importance of this octahedral tilting arrangement we are able to explain why 
the Ca atoms align along (110) planes as opposed to a rock-salt type configuration. Ob- 
servation of the schematics I/II in fig. H] indicates that there can be only one such tilting 
configuration per layer. These configurations must be staggered diagonally in order not to 
impose large distortions on the oxygen octahedra. Such an arrangement excludes the rock 
salt configuration. It also excludes alternating planes composed of Ca and Sr cations. 
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0.756 1.053 0.020 


B 


1.442 


1.025 


1.034 1.067 0.087 


0.646 1.030 0.164 


0.749 1.056 0.030 


C 


1.352 


1.025 


1.058 1.061 0.035 


0.647 1.030 0.162 


0.633 1.052 0.070 



TABLE II. Tabulated data for the A=llollxyl2ol2xy, B=llollxyl2xl2y and C=llollxllyllxy con- 
figurations at 50% Ca, 50% Sr pictured in fig. [TJ The first two columns are the difference in energy 
between the given configuration and the Hollxyl2ol2xy ground state (relative energy) and the c/a 
ratio for the configurations with no rotations or tilting and no FE displacements of the atoms in 
the unit cell. Each of the subsequent sets of three columns give the relative energy, c/a ratio and 
cell polarization with only FE displacements, AFD rotations and tilting about the a-axis of the 
cell. 

In table ITLT1 we compare the lattice parameters and polarizations for the Si(001) lattice 
matched CaTi03 and SrTiOs with the lowest energy configuration of Cao.sSro.sTiOs. The 
c/a ratios sandwich the Cao.5Sr .5Ti0 3 value of 1.043 with 1.002 and 1.053 respectively and 
with respect to the 40 atom cell. We find that the mixed configuration has the lowest net 
polarization. Both Ca- and SrTi03 are considered incipient ferroelectrics in their unstrained 
states but have been found to be polar in their strained configurations^ 1 ^. The oxygen cage 
rotations may be classified as a + b~c - and a°a°c _ for the CaTiOs and SrTiOs respectively. 
In these pure ABO3 perovskites the A-site cations have smaller FE displacements than in 
the mixed case. 

In summary we have identified a ground state configuration for the Si(001) lattice matched 
Ca 05 Sr a5 TiO3. This consists of alternating planes of Ca and Sr running along the (110) 
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Configuration 


a (A) 


b(A) 


c(A) 


Polarization [C/m 2 ] 


Cao.5Sr .5Ti0 3 


7.728 


7.728 


8.06 


0.077 


CaTiC-3 


7.728 


7.728 


7.740 


0.146 


SrTi0 3 


7.728 


7.728 


8.136 


0.304 



TABLE III. Lattice vectors and polarization for Si lattice matched Cao.sSro.sTiOa (lowest energy 
configuration), CaTiOs and SrTi03 as determined in this work. 

direction relative to the cubic ABO3 perovskite phase. We find this structure to be ferroelec- 
tric with a net polarization of 0.077 C/m 2 and a reduced surface strain compared to STO. 
We have identified eight higher energy configurations which, with one exception, exhibit net 
ferroelectric polarizations. The A-site cations counter the dominant ferroelectric contribu- 
tions of the Ti atoms. This work shows that the small ionic radius of the Ca-ions combined 
with octahedral rotations allows them to displace in such a way as to minimize the net cell 
polarization. The placement of the A-site cations determines whether AFD rotations or 
tilting effects are dominant. The former suppresses octahedral distortions and subsequent 
displacement of the A-site cations. When tilting is dominant, the tetragonal distortion is 
enhanced and the cell polarization decreases. 

The present study does not include surface and interface effects which very likely play 
an important role in both the overall structure and polarization of thin film structures. 
However, it identifies the key interactions which govern the stability and polarization of 
AA'TiOs perovskites which should prove valuable in designing such materials. 

Computations were performed at the AFRL DoD Major Shared Resource Center. 
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